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Abstract 

If a black hole starts in a pure quantum state and evaporates completely 
by a unitary process, the von Neumann entropy of the Hawking radiation 
initially increases and then decreases back to zero when the black hole has 
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^f) • disappeared. Here numerical results are given for an approximation to the 
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time dependence of the radiation entropy under an assumption of fast scram- 
bling, for large nonrotating black holes that emit essentially only photons and 
I gravitons. The maximum of the von Neumann entropy then occurs after about 

53.81% of the evaporation time, when the black hole has lost about 40.25% 
of its original Bekenstein-Hawking (BH) entropy (an upper bound for its von 
Neumann entropy) and then has a BH entropy that equals the entropy in the 
radiation, which is about 59.75% of the original BH entropy AitMq, or about 
7.509Mo2 6.268 x W'^'^ {Mq / Mq)^ , using my 1976 calculations that the pho- 
ton and graviton emission process into empty space gives about 1.4847 times 
the BH entropy loss of the black hole. 
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Interest in black hole information [Il|2l[3lll[5l|6l[71[8l[g[T0l[IIl[T2l[l3l[Tl[l5l 
[T6l[I71[T8l[l9l[20l[2ll|2a|2a[2a|25l|2^ 

[MisniiiiiiiiiisiiiaiisiiisiiiziiisiiiniEffl 

[MllMllMlEslEniiZlEHliaiTolini^ 

[85llMl871[M[M[9miOTll92ll9mM[95l[9m[971M 

[insiiinaiinHiiinainniiiiiiiiiainsiinii^ 

[MITiTlfT^ITimiTMfTiSimiTiTimfTiQllCT 

11571 1158[ 11591 116U] has surged recently with the publication by Almheiri, Marolf, 
Polchinski, and Sully |161] of a provocative argument that suggests that an "infalling 
observer burns up at the horizon" of a sufficiently old black hole, so that the horizon 
becomes what they called a "firewall." This paper has elicited a large number of 
responses, some of which support the firewall idea |162[ I163[ I164[ I165[ 1166] . others 
seem rather agnostic |167[ I168[ 1169] , and yet others of which raise skepticism about 

it [1701 [ml [ml [ml [IH [Ea insi 

A central question concerns the entanglement of a black hole with the Hawking 
radiation |184[ 1185] that has been emitted. In fact, it is the argument that at late 
times the Hawking radiation is maximally entangled with what is inside the black 
hole that suggests that what is just inside cannot be significantly entangled with 
what is just outside, and without this latter entanglement, an observer falling into 
the black hole should encounter high-energy radiation at the horizon that would 
burn up the observer |161j . 

It would be impossible for an observer to remain outside the black hole horizon 
and yet to measure directly the entanglement or lack thereof across it. However, if 
the black hole emission is a unitary process, the entanglement between the black 
hole and the earlier Hawking radiation should eventually be transferred to an entan- 
glement between the earlier Hawking radiation and later Hawking radiation. This 
should have an effect on the time dependence of the von Neumann entropy of the 
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radiation up to some point in retarded time during the emission. That would be a 
function of that time that in principle should be calculable (though not yet in prac- 
tice) from a theory of quantum gravity that considers only the black hole exterior 
and indeed only what gets radiated to future null infinity. 

The importance of this question of the time dependence of the Hawking radiation 
entropy has recently been emphasized by Strominger [116] . He gives a brief outline 
of five candidate answers to the question (bad question, information destruction, 
long-lived remnant, non-local remnant, and maximal information return). Here I 
shall assume without proof that the fifth answer is correct and make an additional 
assumption of nearly maximal entanglement between the black hole and the Hawking 
radiation |186[ [25] (with support from the later conjecture of fast scrambling |187[ 
11881 11891 11901 1191] ). In particular, I shall assume that the von Neumann entropy 
of the radiation is very nearly that given by the semiclassical approximation up to 
the point where that equals the Bekenstein-Hawking entropy of the remaining black 
hole, and thereafter is very nearly the same as the Bekenstein-Hawking entropy 
of the shrinking hole. Then I shall give, from my previous numerical calculations 
|192[ I193[ 1194] . an approximation to the time dependence at future null infinity of 
the von Neumann entropy of the Hawking radiation up to that point from a large 
nonrotating black hole. 

I shall consider only nonrotating uncharged (Schwarzschild) black holes that 
are sufficiently large (and hence at low enough Hawking temperature) that they 
emit essentially only massless particles, which I shall assume are only photons and 
gravitons. I shall assume that all other particles (e.g., neutrinos) are massive, that 
the lightest massive particle (perhaps the electron neutrino) has mass m, and that 
one can neglect the radiation from a black hole into particles of mass m if (using 
Planck units) 87r < m/T = 87rmM, where M is the mass of the black hole. My 
numerical calculations |194] showed that for a neutral particle species of spin-half, 
the power emission for m = is only 1.505 x 10"^ times that for spin-half particles 
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with m = 0. so indeed it is a good approximation to neglect the Hawking emission 
of particles with Mm > 1. 

Therefore, I shall restrict attention to large black holes, with M > ^ 
1.336 337 707 x IQ-^^ M^eV / (mc^) [I95]. Any black hole of a solar mass Mq or 
greater would qualify so long as the lightest massive particle has a rest mass energy 
greater than about 0.134 nano-electron volts, which is over ten orders of magnitude 
smaller than the current upper limit of 2 eV on the anti-electron-neutrino mass |195] 
and nearly eight orders of magnitude smaller than 0.0087 eV, the square root of the 
smallest known difference [(7.50 ± 0.20) x 10^^ eV^] of the squares of two of the 
three neutrino mass-energies [195]. None of these data exclude the possibility of 
one massless neutrino species (or of other massless particles besides photons and 
gravitons, provided their coupling to known particles is sufficiently weak), but with- 
out any theoretical reason to expect such massless species, it seems most likely that 
there are no massless particles other than photons and gravitons (for which there 
are strong theoretical reasons from the assumption of gauge invariance for suspect- 
ing that they are precisely massless), and no massive particles of smaller rest-mass 
energy than 1.34 x 10^^° eV. If this is indeed true, as I shall assume here, black holes 
of solar masses and above would radiate almost entirely photons and gravitons. 

My numerical calculations (with an estimated accuracy of about four or five 
digits; not all of the digits given below may be significant, but I give them so that 
the numerical error of my calculations can be found if anyone repeats the calculations 
to higher precision) of the Hawking particle emission rates from both nonrotating 
and rotating black holes |l92l [1931 [M] (see also [1961 ttHl [HH]) in Planck units 
showed (e.g., in Table II of |193j . summarized in |197j and in Eq. (19) of |198j ) that 
the emission power in photons is approximately 0.000 033 638/M^ and in gravitons 
is close to 0.000 003 836/M^, for a total emission rate of massless radiation energy 
of a/M^ ^ 0.000 037474/M2 with 

a = -M^^ ^ 0.000 037474. (1) 
dt ^ ' 
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My numerical calculations also showed (e.g., summarized in |197] and in Eq. (20) 
of |198j ) that the rate of entropy emission in photons is approximately 0.001 268 4/M 
and in gravitons is about 0.000 130 0/M, for a total emission rate of massless radi- 
ation entropy of dSj-a^d/dt ~ 0.001 398 4/M. These latter numbers were calculated 
from the former by using the results I reported in the last paragraph of the first 
column of the tenth page (page 3269) of [193]: 

"There is still some entropy produced by the partial scattering off the gravita- 
tional potential barrier surrounding the hole, but outside the supperradiant regime 
this can only partially cancel the entropy flow out of the hole and serves in effect 
to increase the entropy emitted to the surrounding region for a given entropy loss 
by the hole. For example, numerical calculations for a nonrotating hole show that 
the emission of s = | particles into empty space increases the external entropy by 
1.6391 times the entropy drawn out of the hole, s = 1 particles increase it by a 
factor of 1.5003, s = 2 particles by 1.3481, and the canonical combination of species 
[the set of species known in 1976 with masses less than 20 MeV: gravitons, photons, 
electrons and muon neutrinos with one helicity each, electrons, and the correspond- 
ing antileptons, all taken to have Mm <^ 1 for this numerical result] gives 1.6233 
times as much entropy in radiation as the entropy decrease of the hole." 

When one uses the fact that the Bekenstein-Hawking semiclassical entropy |199[ 
I200[ 1185] of a nonrotating uncharged black hole is, in Planck units, ^bh = 47rM^ 
|185] . the rate of entropy decrease of the hole is Sna/M. From just photon emission, 
one would get (cf. [HT] and Eq. (20) of [ISH]) 0.000 845 41/M, and from just graviton 
emission one would get 0.000 096 41/M, for a total rate of decrease of the Bekenstein- 
Hawking black hole entropy of -dS^Yi/dt = Sica/M ^ 0.000 941 82/M. This then 
gives the ratio by which the coarse-grained entropy of the Hawking radiation (e.g., 
ignoring its entanglement with the black hole) is greater than the coarse-grained 
Bekenstein-Hawking entropy of the black hole as 

^f-'^/f^^ ^ 1.484 72; (2) 
—dSBu/dt 



the last one or two digits probably are not significant but are given as a challenge 
for someone to do a higher-precision calculation to find the error. (This constant P 
should not be confused with the (3 given in Table II of |193] . which gave a measure 
of the angular momentum emission of very slowly rotating black holes.) 

From the fact that the mass of a large nonrotating black hole that emits al- 
most entirely just massless photons and gravitons decreases according to dM/dt = 
—a/M^, one can get that the time evolution of the black hole mass is 

M{t) = {Ml - Saty/' = [3a(tdecay - t)]'/' = Mo(l - t/tdccay)'/', (3) 

where Mq is the initial mass of the black hole at time t = and where 

1 / M \ ^ 

tdecay = jM^ = —M^ ^ SSOSM^ ^ 1.159 X 10^^ f ^ j yr (4) 

is what the decay time would be of a nonrotating uncharged black hole of initial 
mass Mo if only massless photons and gravitons were emitted, with 

^ " ^ " -3MW/dt " 
Assuming that the lightest massive particle has mass m, the emission of massive 
particles would reduce the actual decay time by an amount of the order of 1/m^ in 
Planck units, e.g., by a time of the order of 10^^(eV/m)^ years, which for m ~ 0.1 
GeV is greater than the age of the universe but which for m ^ 10^^° eV (which I 
am assuming) is much less than the decay time tdccay ~ 10^'' years for a solar-mass 
black hole, so in this paper I shall ignore any such correction to the black hole decay 
time. Another correction that I shall mostly ignore but shall return to later is the 
uncertainty of the decay time from the stochastic nature of the Hawking emission of 
very roughly ~ Mq particles, which would lead to fluctuations in the total decay 
time of an ensemble of black holes of initial mass Mq that would be expected to be 
of the order of 6t ~ tdccay/ViV ~ tdccay/Mo ~ M^ ~ lO'^^iMo/M^y years. These 
fluctuations in the total decay time have magnitudes that are comparable (with 
c = 1) to the fluctuations in the position of the final decay of the black hole in the 
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center-of-mass frame of the initial black hole, 5x ~ Mq ~ 1Q'^^{Mq/ MqY light-years 
[21 [1571 [TS]. 

The semiclassical approximation for the black hole emission |184[ I185[ [T] would 
then give the time- dependence of the entropy of a decaying large nonrotating un- 
charged black hole and of the emitted radiation as 

/ t ( t \'/' 

^BH(t) = 4vrMM 1 - — 3 -47rMoMl-— — 3 , (6) 



7M^ J ° V 8895MoV 
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} 7M0V 



t \ 



2/3 



mM§J 



^ 47r(1.4847)M2 

(7) 

Here I have put tildes on the expressions for these entropies to emphasize that 
they are the entropies obtained from semiclassical approximations. Ultimately we 
would like instead to obtain von Neumann entropies of black hole and radiation 
subsystems (S'vn(^)) of the entire universe, which for simplicity I shall assume is in 
a pure quantum state with zero von Neumann entropy. Also for simplicity I shall 
assume that one is considering a black hole that has formed in a pure quantum state, 
which implies that it has no quantum entanglement with the rest of the universe. 
Of course, these are unrealistic assumptions, but they might give a reasonably good 
approximation if one focuses on one pure component of the quantum state of the 
universe in which the black hole forms from the gravitational collapse of a star with 
a von Neumann entropy of the order of the number of particles in the star, say 
S ~ 10^^, which is much smaller than the Bekenstein-Hawking entropy of a stellar 
mass black hole that is at least S'bh = 47rMQ ~ 10''^. In other words, in comparison 
with the von Neumann entropies to be discussed that are at least of the order of 
10''^, I shall neglect the much smaller ordinary stellar entropies that are of the order 
of 10^^. 

I shall take the semiclassical entropies to be upper limits on the von Neumann 
entropies of the corresponding subsystems with the same macroscopic parameters. 
That is, I shall first assume that the semiclassical Bekenstein-Hawking entropy of 
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a nonrotating uncharged black hole is a good approximation for the maximum von 
Neumann entropy of a black hole of the same energy, at least if one neglects entropy 
associated with the location and/or motion of the black hole in a space sufficiently 
large that this could in principle rival the Bekenstein-Hawking entropy. Secondly, I 
shall assume that the semiclassical entropy calculated for the Hawking radiation is 
the maximum von Neumann entropy for radiation with the same expectation val- 
ues of the numbers of particles in each of the modes. This second assumption is 
supported by the fact that indeed the semiclassical approximation does give a ther- 
mal density matrix for each radiation mode, and no entanglement between different 
modes, though because of the black-hole greybody factors the temperature varies 
from mode to mode. 

Now under the assumptions that the black hole starts in essentially a pure state 
(von Neumann entropy much less than the Bekenstein-Hawking entropy), that the 
Hawking evaporation is a unitary process, and that we can neglect the interaction 
with other systems (e.g., no large external source of entropy impinging upon the 
black hole or interacting with the outgoing radiation), the von Neumann entropy 
of the evaporating black hole equals that of the Hawking radiation that has been 
emitted. When combined with the assumptions, of the previous paragraph, we reach 
the conclusion that at retarded time t, the von Neumann entropy of the black hole 
or of the Hawking radiation cannot exceed the minimum of the two semiclassical 
entropies given above, 5'bh(^) and 5'rad(^)- 

Next, based upon what is suggested by the results of [1861 EHl [Ml IHHl EM IM 
1191] ). as discussed above, I shall assume that the actual von Neumann entropy of 
the Hawking radiation up to retarded time t (which is also that of the remaining 
black hole, though from now on I shall focus on the radiation entropy as potentially 
being more easily measured) is very near to its maximal allowed value that is the 
minimum of Sbh(^) and S'rad(^)- 

'S'bh(^) decreases monotonically with time for < t < tdecay = 7^o ~ 8895Mq, 
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and S'rad(i) increases monotonically with time. The two values cross at the time 

3/2' 



1 



■^decay — ^^decay ~ 0.53810 t^ecay 



4786M3 



with 
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at which time the mass of the black hole is 

1/2 



4 1 

/3 



0.53810, 



/3 + 1 



O.773OIM0, 



(8) 



(9) 



(10) 



and the entropy of the radiation and of the black hole is 



5** = 5'vN(i*) — 'S'BH(i^*) — 'S'i.ad('^*) 



/3 



^^-j 47rMo^^ 0.59754 5bh(0) 
= 0.59754(47rMo2) ^ 7.5089M2 ^ 6.268 x lO^^(Mo/M0)l (11) 

This time, t* ~ 0.5381 tdecay, under the assumptions being made in this paper, 
would be the retarded time at which the von Neumann entropy of the Hawking 
radiation from a large nonrotating uncharged black hole, which was initially in 
essentially a pure quantum state, reaches its peak value and thereafter decreases. 
It is what I think is intended to be meant by what is called the 'Page time' in the 
literature [liIl[l6ll[T62l[I70l[I72l[T63l[T6^ 

fT59j . though often [Ml [1621 [TTOl [1631 [M [1571 [M] it is somewhat inaccurately 
described as something equivalent to the time "when the black hole has emitted 
half of its initial Bekenstein-Hawking entropy" |161] . which would be correct if /3, 
the ratio of the semiclassical entropy emitted in Hawking radiation to the decrease 
in the Bekenstein-Hawking semiclassical entropy of the black hole, were unity, but 
it is somewhat different when /3 7^ 1, as for the emission of photons and gravitons 
from a nonrotating uncharged black hole that gives /3 ~ 1.4847. In particular, 
the time at which the black hole area has decreased to half its original value is 



t 



1/2 



1 - 2-3/2)t, 



decay 



0.64641 



decay 



1.201 1*. 



As we can sec above, under the assumptions that lead to the von Neumann 
entropy of the radiation being (very nearly) the minimum of the semiclassical en- 
tropies of the black hole and of the Hawking radiation, the time is actually about 
53.81% of the black hole lifetime (or about 83.24% of the time until half the area 
and half the Bckcnstcin-Hawking entropy are lost), but on the other hand at this 
time the black hole has lost only about 40.25% of its original Bekenstein-Hawking 
entropy. The excess entropy in the radiation over what is lost from the black hole 
comes from the semiclassical entropy generated in the emission process, which is 
a nonequilibrium transfer of energy from a black hole of positive temperature to 
originally empty space of zero temperature. That is, the loss of about 40.25% of 
the original black hole entropy generates radiation of /3 ~ 1.4847 times as much, or 
about 59.75% of the original black hole entropy, equahng the remaining black hole 
Bekenstein-Hawking entropy. 

Now under the assumptions above, we can say that the von Neumann entropy of 
the Hawking radiation Sy^{t), as a function of the retarded time t, is very nearly the 
semiclassical radiation entropy S'rad(^) for t < and is very nearly the Bekenstein- 
Hawking semiclassical black hole entropy ^bh(^) for t > t^. Using the Heaviside 
step function 9{x), the von Neumann entropy of the Hawking radiation from a large 
nonrotating uncharged black hole may be written as 



-(- 



t 



) 



2/3" 




Sy^{t) ^ 4:71/3 1 



e{u -t) + All Ml 1 



1 - - — e{t - u) 



lecay 



^decay 




2/3" 



fa 47r(1.4847)Mo2 1-1 



8895Mo^ 

2/3 



^(4786M3 - 1) 



e{t- 4786Mo^) 



(12) 



See Figure 1 for a plot of this entropy versus retarded time. 
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Figure 1: Hawking Radiation Entropy vs. Time 
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Because of the fluctuations in the black hole emission rate that make the time 
for a black hole to evaporate down from an initial mass of Mq to M^, ^ 0.7730Mo 
uncertain by an amount of the order of Mq in Planck units, one would expect the 
peak in the von Neumann entropy function S^sit) at t = to be rounded off over 
a time scale of the order of Mq, which of course for large black holes is very much 
smaller than ^ 4786Mq, smaller by a factor of the order oH/Mq that is ~ 10^^^ for 
a solar mass black hole. Furthermore, if the von Neumann entropy of the radiation 
up to retarded time t is evaluated at future null infinity up to its intersection with a 
null cone moving outward from the center of mass position of the total system of the 
black hole plus radiation in the center-of-momentum frame of that entire system, 
the fluctuations of the black hole position from that center of mass that are of the 
order of Mq [21 11571 1181] would also be expected to contribute to a rounding off of 
the entropy function S'vn(^) by a comparable amount. Since \dSv^{t)/dt\ ~ 1/Mq on 
either side of the idealized sharp peak (ignoring this rounding off), the rounding off 
over a time of the order of Mq would be expected to reduce the maximum value of the 
peak by an amount ~ Mq. Therefore, I would expect the peak of the von Neumann 
entropy of the Hawking radiation actually to be 5** ^ [47r/3/(/3 + l)]Mo — 0(Mo) ^ 
0.59754 ^bh(O) = 0.59 754 (47rM2) ^ 7.5089M2 ^ 6.268 x W^Mq/Mq)^, with a 
correction (beyond the numerical errors of the coefficients) only of the order of 

Mq ~ W^\Mo/Mq). 

One might think that the fluctuations of the black hole position would greatly 
increase the von Neumann entropy of the Hawking radiation, since if one uses ra- 
diation wavepacket modes of definite angular momenta relative to the center of 
mass of the entire system, the deviation of the black hole from the center of mass 
when it emits a mode will cause it to have a typical angular momentum of the or- 
der of the frequency of the mode, ~ 1/M ~ I/Mq, multiplied by the deviation of 
the black hole position from the center of mass, ~ Mq, for an angular momentum 
j ~ Mq about the center of mass. For each time period ~ M ~ Mq over which 
roughly one particle is emitted by the black hole, the number of angular momen- 
tum modes into which the particle could be emitted is ~ ~ Mq. Therefore, the 
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total emission of ~ Mg Hawking particles would be spread over ~ Mq modes quasi- 
localized in both time and angular momentum, with each mode having a probability 
of being occupied that is p ~ M^"^ and hence a von Neumann entropy per mode 
'S'mode ~ —p Inp ~ Mq'^ In Mq. If one then sums this von Neumann entropy per mode 
over the ~ Mq modes that each have roughly this entropy, one would get a total 
entropy in the radiation of S^ad with flucts ~ •S'mode ~ Mq In Mq, roughly a factor of 
the order of IuMq ~ 87 + ItlIMq/Mq) times the Bekenstein-Hawking semiclassical 
entropy of the original black hole of mass Mq. 

However, this entropy calculated for the Hawking radiation coming from a black 
hole, with fluctuations in its position ~ Mq from fluctuations in the momenta of the 
previously emitted quanta, is a highly coarse-grained entropy, ignoring the quantum 
entanglement between the different modes. For example, the apparent source loca- 
tion for a particle emitted at late times will be roughly where the black hole was at 
the emission time, which depends on the history of the emission of momenta of the 
earlier particles that determines the black hole recoil that leads to the fluctuations 
in position. When this quantum entanglement between modes is taken into account 
(which is absent in the semiclassical approximation that has the black hole remain 
at the flxed location of the center of mass of the total system, that of the original 
black hole before any radiation was emitted, and decay with a fixed M{t) that just 
depends on the expectation value of the emission rate with no quantum fluctuations 
of the black hole in the semiclassical approximation), I would expect that the von 
Neumann entropy of the radiation is close to that of the semiclassical approxima- 
tion up to the time at which the semiclassical entropy of the radiation reaches 
the Bekenstein-Hawking entropy of the black hole, which is presumably a good ap- 
proximation for the maximum von Neumann entropy of the black hole that always 
stays equal to the von Neumann entropy of the radiation under my assumptions of 
a pure initial state, unitary evolution, and the possibility to divide the total system 
up into black hole and radiation subsystems. 

Now the stochastic recoil of the black hole might be expected to increase its von 
Neumann entropy (and hence also that of the Hawking radiation), but this appears 
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to be by an amount that is quite negligible compared with the Bekenstein-Hawking 
entropy ~ Mq. In particular, since the black hole develops a position uncertainty ~ 
Mq from its momentum uncertainty ~ 1 (all in Planck units) from the square root of 

~ Mq particles emitted each with momentum ~ I/Mq, the number of momentum 
states up to momentum ~ 1 in the region of volume ~ Mq would be of the order of 
~ Mq , which would lead to an entropy of the order of the logarithm of this number 
of states, or ~ 61nMo ~ 522 + 61n(Mo/M0). Although this is a large compared 
with unity, and would lead to a number of states far larger than a googol raised to 
the nine-fourths power, it is miniscule compared with the maximum von Neumann 
entropy of the radiation that is here estimated to be 5* ~ 6.268 x 10'^^{Mq/Mq)'^, 
and also very much smaller than the uncertainty of this that I have estimated to be 
~ Mo ~ 1038(Mo/M,). 

In conclusion, under the assumptions that a large (say solar mass or greater) 
nonrotating uncharged (Schwarzschild) black hole starts in nearly a pure quan- 
tum state and decays away completely by a unitary process while being nearly 
maximally scrambled at all times, the von Neumann entropy of the Hawking ra- 
diation (almost entirely photons and gravitons for such a large black hole, as- 
suming no other particles of rest mass m ~ 10^^° eV) increases up to a max- 
imum of ^ 0.59754(47rM2) ^ 7.5089M2 ^ 6.268 x lO^'^(Mo/M0)2 at time 
t, ^ 0.53810 tdecay ~ 4786Mq3 ^ 6.236 x 10^'^{Mo/MQ fjr and then decreases back 
down to near zero (the von Neumann entropy of the original black hole, assumed to 
be small). The numerically estimated time dependence of the von Neumann entropy 
of the Hawking radiation is given by Eq. ( TT2l) and is plotted in Figure 1. 
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ber 30 - December 1 firewall conference at Stanford University, and to the partic- 
ipants there (particular Andrew Strominger, who asked explicitly about the time 
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entropy of the Hawking radiation. This work was supported in part by the Natural 
Sciences and Engineering Research Council of Canada. 
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